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ABSTRACT 

Aims. For studying how the field lines are twisting nearby the light cylinder surface, which provides the free motion of AGN winds 

through the mentioned area, the investigation of the centrifugally driven curvature drift instability is presented. 

Methods. Studying the dynamics of the relativistic MHD flow close to the light cylinder surface, by applying a single particle approach 

based on the centrifugal acceleration, the dispersion relation of the instability is derived and analytically solved. 

Results. Considering the typical values of AGN winds, it is shown that the time scale of the curvature drift instability is much less 

than the accretion process time scale, indicating that the present instability is very efficient and might strongly influence processes in 

AGN plasmas. 
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1. Introduction 

For studying the AGN winds the fundamental problem relates to 
the understanding of a question: how the plasma goes through 
the Light Cylinder Surface (LCS - the hypothetical zone, where 
the linear velocity of rotation equals the speed of light). An in- 
nermost region of AGNs is characterized by the rotational mo- 
tion, and it is obvious that such a character of motion must affect 
the plasma dynamics. According to the standard model, the mag- 
netic field due to the frozen-in condition undergoes plasmas and 
consequently the flow follows the field lines, co-rotating with 
them. This means that the plasma particles moving along quasi 
straight magnetic field lines in the nearby area of the LCS have 
to reach the speed of light. On the other hand no physical sys- 
tem can maintain such a motion and this fact must result in a 
certain twisting process of the magnetic field lines on the LCS. 
Generalizing the work: (Machabeli & Rogava 1994) for curved 
trajectories, by Rogava et al. (2003) the dynamics of a single 
particle moving along a prescribed rotating curved channel has 
been investigated. It was shown that if the trajectories are given 
by the Archimedes spiral, the particles can cross the LCS avoid- 
ing the light cylinder problem. Making one more step to this 
investigation is to find an appropriate mechanism, which might 
provide the twisting of the magnetic field lines, giving rise to the 
shape of the Archimedes spiral, which in turn makes the dynam- 
ics force-free. 

The light cylinder problem in the context of the force- 
free regime has been extensively studied numerically for pul- 
sars. The investigation developed in (Spitkovsky & Arons 2002, 



Spitk ovsky 2004) has shown that the plasmas can go through the 



LCS. This work was based on a current generated by the elec- 
tric drift (Blandford 2002), which is vanishing for quasi neutral 
plasmas and can not give a contribution in the dynamics of as- 
trophysical flows having almost equal numbers of positive and 
negative charges. 



Since the innermost region of AGNs is rotating, a role of 
the Centrifugal Force (CF) seems to be particularly interesting 
for the study of the relativistic plasma motion. The centrifugally 
driven outflows have been discussed in series of works. A spe- 
cial attention deserves the work of Blandf ord & Payne (1982)| 
where the authors discuss the possibility of the energy and 
angular momentum pumping from the accretion disk, empha- 
sizing the role of the centrifugal acceleration in this pro- 
cess. By [Gangadha ra & Lesch 1997| the CF was considered in 
the context of the non thermal radiation from the spinning 
AGNs. Generalizing the mentioned work it has been shown 
dOsmanov et al. 20071|Rieger & Aharonian 2008| > that due to the 
centrifugal acceleration, electrons gain very high energies with 
Lorentz factors up to y ~ 10 8 . This means that the energy budget 
in the AGN winds is very high and if one finds a mechanism of 
the conversion of at least a small fraction of this energy into a 
variety instabilities, one might have interesting consequences in 
the physics of AGN outflows. 

The centrifugal force may drive different kinds of instabili- 
ties. Obviously the CF acting on a moving particle changes in 
time and in the context of instabilities plays a role of the pa- 
rameter. Consequently the corresponding instability is called the 
parametric instability. 

The centrifugally driven parametric instability first has been 
introduced in ( Machabe li et al. 20051 1 for the Crab pulsar mag- 
netosphere. We have argued that the centrifugal force may cause 
the separation of charges, leading to the creation of an unsta- 
ble electrostatic field. Estimating the linear growth rate it has 
been shown that the instability was extremely efficient. The 
method developed in ([Machabeli et al. 20051 ) was applied for 
AGN jets dOsmanov 20081 1 for studying the stability problem 
of the rotation induced electrostatic instability and for under- 
standing how efficient is the centrifugal acceleration in this pro- 
cess. Another kind of the instability which might be induced by 
the CF is the so called Curvature Drift Instability (CDI). Even 
if the field lines initially have a very small curvature, it might 
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a - {e, b], 



where 




Fig. 1. Two orthonormal bases are considered: i) cylindrical 
components of unit vectors, (e$, e^, e A ); ii) unit vectors of the 
system rigidly fixed on each point of the curve, (e r , eg, ej, re- 
spectively. C is the center of the curvature. Hereafter the set of 
coordinates I call the field line coordinates. 



cause a drifting process of plasmas, leading to the CDI. By 
Osmanov et al. (2008a) the two component relativistic plasma 
has been considered for studying the role of the centrifugal ac- 
celeration in the curvature drift instability for pulsar magneto- 
spheres. The investigation has shown that the growth rate was 
more than pulsar spin down rates by many orders of magnitude, 
indicating high efficiency of the CDI. The curvature drift current 
produces the toroidal component of the magnetic field, which 
due to the efficient unstable character of the process amplifies 
rapidly, changing the overall configuration of the magnetic field. 
This leads to the transformation of field lines into the shape of 
the Archimedes spiral, when the motion of the particles switches 
to the so called force-free regime (Osmanov et al. 2008b) and 
the plasma goes through the LCS. 

In the present paper in order to investigate the twisting pro- 
cess of magnetic field lines due to the CDI, the method devel- 
oped in (Osmanov et al. 2008a, Osm anov et al. 2008bl > will be 
implemented for AGN winds. 

The paper is arranged as follows. In ^2]we introduce the cur- 
vature drift waves and derive the dispersion relation. In |3] the 
results for typical AGNs are present and in $4] we summarize 
our results. 



2. Main consideration 

We start the investigation by considering the two com- 
ponent plasma consisting of the relativistic electrons with 
the Lorentz factor y e ~ 10 5 ~ 8 dOsmanov et al. 2007 , 
Rieger & Aharonia n 2008[ ) and the bulk component (protons) 
with jb ~ 10. Since we are interested in the twisting process, 
we suppose that initially the field lines are almost rectilinear in 
order to study how this configuration changes in time. 

Let us start by the Euler equation governing the dynamics 
of plasma particles, co-rotating with the straight magnetic field 
lines. By applying the method developed by Che dia et al. (1 996) 
one can show that the Euler equation gets the following form: 



+ (v ff V)p„ 

at 



- r ^+^(E + v„xB), 

m„ 



(1) 



Vl - n 2 R 2 /c 2 . 



Here p ff is the momentum, y a - the velocity and y a - the 
Lorentz factor of the relativistic particles. E and B are the elec- 
tric field and the magnetic induction respectively. By q a and m a 
we denote particle's charge and the rest mass respectively. We 
express the equation of motion in the cylindrical coordinates (see 
Fig. [TJ. The first term of the right hand side of the Euler equa- 
tion — y ff £V£ represents the centrifugal force, which on the light 
cylinder surface becomes infinity. This means that its overall ef- 
fect is significant in the nearby zone of the LCS. For describing 
our physical system, one needs the full set of equations, and in 
order to close the system, we add to Eq. (Q]) the continuity equa- 
tion: 



dn <* it/ \ n 

— + V(n a \ a ) = 0, 
at 



and the induction equation: 
1 <9E 

c dt c 



„ „ 1 8E 4-7T V~! 

a-e,b 



(2) 



(3) 



where n a and J„ are the density and the current, respectively. 

In the zeroth approximation the plasma particles undergo 
only the centrifugal force. Different species at different posi- 
tions experience different CFs, which will cause the separation 
of charges, leading to the creation of the additional electromag- 
netic field considered as the first order term in our equations. 

The leading state is characterized by the frozen-in condition, 
E + v a o x Bo = 0, which reduces Eq. ([TJ into the following form 
(|Machabeli & Rogava 1994J: 



dv 
dt 



D. 2 R 



1 



n 2 R 2 



2p2 



1 - 



&R 



2y2 

A 



(4) 



The present equation of motion describing the kinematic be- 
haviour of the single co-rotating particle has the following solu- 
tion: 



v(f) = v., ~ c cos(Qf), 



(5) 



for ultra relativistic cases (y » 1) and the following initial con- 
ditions: R(Q) - 0, v(0) w c (v = dR/dt). Here v„ denotes the 
velocity component along the magnetic field lines. 

Since we suppose that the magnetic field lines initially have 
the small curvature, the particles moving radially, will drift along 
the x axis as well (Osmanov et al. 2008a) (see Fig.[TJ. The men- 
tioned drift of charges will produce the corresponding current, 
which inevitably will create the toroidal magnetic field, chang- 
ing the overall configuration of the field lines. Therefore the aim 
of the present work is to study the role of the centrifugaliy in- 
duced curvature drift instability in the twisting process of the 
magnetic field lines. For this purpose one can linearize the sys- 
tem of equations Eqs. (Q]|3]l, perturbing all physical quantities 
around the leading state: 



T : 



,v,p,E,B) 



(6) 
(7) 
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Let us express the perturbation by following: 

x i' l (t,r)cc y ¥ l (t)exp[i(kr)] , (8) 

then from Eqs. (Q]|3]l one can derive the linearized set of equa- 
tions governing the CDI: 



- i(k x u a + k^v^Pax = — v„Bj, 

Of f^ffr 



dn l 

- i(k x u a + k^v t )n a = ik x n a v ax , 
- ik^cB]. = An ^ q a (n° a v l ax + n l a u a ). 



(9) 

(10) 
(11) 



Here by u a we denote the curvature drift velocity along the x 
axis: 



Xln — 



y°» v i 

o>bRb 



(12) 



where <y [vB = q a B {) /m a c; R B is the curvature radius of magnetic 
field lines; and 



Bo 



2L 



is the equipartition magnetic induction on the LCS for the lead- 
ing state (here L is the luminosity of the AGN and R/ c — c/Q 



is the light cylinder radius). For deriving Eqs. (I9¥1T|> the wave 
propagating almost perpendicular to the equatorial plane has 
been considered and the expression: v' * cE\/Bq was taken into 
account. For simplicity the set of equations are given in the co- 
ordinates of the field line (see Fig.Q]). 

Let us express v ax and n\ by the following way: 



v« = V ax e ti 



where 



A ax (t) = -j- + ^ sin(20f), 



A a<p (t) = — sin(Qf). 



(14) 
(15) 

(16) 
(17) 



Then, by substituting Eqs. (1 141 15b into Eqs. (l9lfTT1i. it is easy to 
solve the system for the toroidal component: 



-ik.cBlit) = V ^e ikA «« f 

I y fi^e'WO f ' dt > f 



e -*A„(O v < t ')B r {t')dt'+ 



-^V\(t")B r (t")dt" 



(18) 



where w„ = e ^Ann^lm a is the plasma frequency. After making 
the Fourier transform (see Appendix A), Eq. ( fT8l gets the form: 



B 



J s (ga)J„(h)Ji(g a )J p (h) y 



xB r (oi + Q (2[s - /] + n - p + a)) 

to 2 a k x u a 



1 



/l %f£(y 



Q(2s + n) 

JAga)Jn{h)Jl{g a )J p {h) 



z z 



a= e ,4 V cr/i=±1 Sj „ jlp (a) + h^L + Q(2[j + ^] + «)) 
x£ r (w + Q(2[s-Z+yu] +n-p + o-)), (19) 



where 



k x ii a 
~4Q' 



h - ^ C 

1 ~ ~o~' 



Let us note that Eq. ( fT9] l is written for B r (a>) and on the 
right hand side of the equation, there is the infinite number of 
components with B,{co + Q), B r (a> + 2Q),... etc. This means that 
for solving the mentioned equation, thus for closing the system, 
one needs to add the corresponding expressions for B,(co ± Q), 
B,-(aj + 2Q.),... etc. But then the system becomes composed of 
the infinite number of equations, making the task unsolvable. In 
order to overcome this problem one has to use a certain, phys- 
ically reasonable cutoff on the right hand side of the equation 
dSilin et al. 19701 I. As it is clear from Eq. ( fT9l , the considered 
instability is characterized by the following proper frequency of 
the curvature drift modes: 



(13) co 



k \f&(y 



(20) 



when the corresponding conditions: k x u a /2 < 0, 2s + n — and 
2[s + fi] + n — are satisfied. As we will see only the resonance 
terms give the significant contribution to the result. 

Let us consider parameters L ~ 10 44 erg/s, Q = 3 x 10 s -1 , 
y e o ~ 10 5 , R B ~ R k ~ 0.001c»r 3 typical for AGN 
winds. Then examining the curvature drift waves with A ~ R[ c 
(A = 2n/k is the wave length), one can show that \k x u e /2\ ~ 
lO -12 ^ 1 << Q (here, it is supposed that k x < and u e o > 0, 
otherwise the frequency becomes negative). Therefore all terms 
with non zero Q.(2s + n) and Q(2[s+fi] +ri) are rapidly oscillative 
and do not contribute in the final result. Consequently the only 
terms, which influence the solution of Eq. ( fT9b are the leading 
terms, the contribution of which reduces the equation (Osmanov 
et al.2008a) (see Appendix B): 



* ' 0~,jU=±l sj 



where 



Hq - 2S ± i - 



tojk x u e 
2y epo k g c' 



(21) 



(22) 



n' = -2(s + p), p'=-2l + LT, g = ge- 

Note, that Eq. (f2TT > does not consist of the bulk components at all, 
because as a direct calculation shows, their contribution is neg- 
ligible with respect to the terms corresponding to the relativistic 
electrons. Expressing the frequency by the real and imaginary 
parts: ai = a>o + iT it is straightforward to show from Eq. (fST), 
that the increment of the CDI writes as : 



2 2j ^ J -^8)J-2( S +M)( h )Ji(g)J-2i+Ah) 



0",//=±l s,l 



(23) 
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Fig. 2. The dependence of logarithm of the instability time scale 
on the normalized wave length. The set of parameters is: y e o = 
{10 5 ; 10 6 ; 10 7 ; 10 8 },% ~ Ri c ,n e0 = 0.001cm- 3 , 1$ = 100% and 
L/L E = 0.01. 
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eO J 



Fig. 3. The dependence of logarithm of the instability time scale 
on the density normalized by the medium density. The set of 
parameters is: y e0 = {10 s ; 10 6 ; 10 7 ; 10 8 }, R B ~ %, A e0 = %, 
= 100% and L/L E = 0.01. 

3. Results 

In this section we investigate the efficiency of the CDI for AGN 
winds. The behaviour of the growth rate will be studied ver- 
sus the wave length, the density of relativistic electrons, their 
Lorentz factors and the AGN bolometric luminosity. 

For studying the behaviour of the instability as a function of 
the wave length one can examine the typical AGN parameters: 
M BH = 10 8 x M Q , Q = 5 x lO^V 1 and L = 10 44 erg/s, where 
Mbh is the AGN mass, M - solar mass and L is the bolometric 
luminosity of the AGN. 

Let us consider Eq. d23l and plot logarithm of the instability 
time scale 

t = I (24) 



10.5 



10 



9.5 



8.5 



8 



7.5 





- 
















^\L/L E =0.01^^\^^ ^^^^ 







5 5.5 6 6.5 7 7.5 8 



Fig. 4. The dependence of logarithm of the instability time scale 
on log{y e o). The set of parameters is: % ~ Rk, n e o/n m - 0.001, 
A e0 = Ric, A^ = 100% and L/L E = {0.01; 0.05; 0.3; 1). 

versus the wave length normalized by the light cylinder radius. 
The present consideration is based on the centrifugal accelera- 
tion. As it has been shown in (Osman ov et al. 2007K due to the 
CF, the relativistic particles may reach very high Lorentz fac- 
tors. For this purpose it is reasonable to investigate the efficiency 
of the instability versus the wave length but for different values 
of Lorentz factors. Fig. [2] shows the mentioned behaviour for 
the following parameters: y M = {10 5 ; 10 6 ; 10 7 ; 10 8 }, R B ~ R lc , 
n M = 0.001cm~ 3 , A^ = 2n/k^ = 100% and L/L E = 0.01, where 
L E = I0 46 erg/s is the Eddington luminosity for the given AGN 
mass. Different curves correspond to different values of Lorentz 
factors. As it is clear from the plots, the time scale is a con- 
tinuously increasing function of A x (= 2n/k x ), which is a direct 
consequence of Eqs. (1221231 1. Indeed, as we see from Eq. (l22l . 
Hrj,±i ~ lMx, which combining with Eq. d23l gives r ~ 1/ y[A~ x . 
Therefore the bigger the initial perturbation wave length, the less 
the instability time scale and consequently the less the CDI effi- 
ciency. For the given range of A x and different values of y e o, the 
CDI time scale varies from ~ 10 7 s {A K /Ri c = 1, y e o = 10 8 ) to 
~ 10 9 (AJRic = 2, r e0 = 10 5 ). 

Sine the CDI growth rate depends on the plasma frequency 
[see Eq. ( 123I22H . which in turn is the function of the density, it 
is obvious that the instability time scale must be influenced by 
the density of relativistic electrons in the AGN winds. Indeed, 
In Fig. [3] the plots of log(f) versus the AGN wind density are 
shown and one can see that the time scale is the continuously 
decreasing function of n e o/n m . The set of parameters is the same 
as in Fig. [2] except A x = R/ c and n e o/n m € {0.001, 1). Here n e o is 
normalized by the intergalactic medium density, n m ~ lcm' 3 . As 
we see from the figure, r varies from ~ 10 9 s (n e o/n m = 0.001, 
y e0 = 10 5 ) to ~ 10 6 s (n e0 /n m = 1, y e0 = 10 8 ). 

In Fig. |4]the plot of log(f) versus log(y e o) is shown for dif- 
ferent luminosities. The set of parameters is the same as in Fig. [2] 
except the continuous range of y e o £ { 10 5 ; 10 8 }, different values 
of the luminosity L/L E = {0.01; 0.05; 0.3; 1} and A x = R, c . The 
figure shows the continuously decreasing behaviour of the insta- 
bility time scale, which is a natural consequence of the fact that 
the more energetic electrons will induce the curvature drift insta- 
bility more efficiently. Indeed, from Eq. (fT2l it is clear that the 
drift velocity is proportional to the Lorentz factor of the particle, 
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Fig. 5. The dependence of logarithm of the instability time scale 
on log(y e Q). The set of parameters is: Rg ~ Ri c , n e o/n„, - 
{0.001; 0.01; 0.1; l},^o = Ric, <L = i0QR ic and L/L E =0.01. 



and hence the corresponding instability will be more efficient, 
leading to the decreasing behaviour of log{r). As we see the in- 
stability time scale varies from ~ 10 10 s (y e o = 10 5 , L/Le = 1) to 
~ 10 8 s (y e o = 10 7 , L/L E = 0.01). On the other hand the plots for 
different luminosities exhibit another feature of the behaviour of 
t: by increasing the luminosity of the AGN, the corresponding 
instability becomes less efficient. 

To see this particular feature more clearly, let us look at Fig. [5] 
exhibiting the dependence of log(f) on L/L E for different values 
of densities. From the plots it is seen that by increasing the lumi- 
nosity, the time scale continuously increases. This behaviour fol- 
lows from the fact, that the bigger the luminosity, the bigger the 
magnetic field [see Eq. ([T3l l and hence the less the drift velocity, 
leading to the less efficient CDI shown in the figure. By consid- 
ering bigger values of densities, the CDI becomes more efficient, 
which we have already explained while considering Fig. [3] For 
the mentioned area of quantities (see Fig.[5j, the time scale varies 
from ~ 10 7 s (L/L E = 0.01, n e0 /nj = 1) to ~ 10 I0 s (L/L E = 1, 
n e0 /nj = 0.001). 

We see from the present investigation that the instability time 
scale varies in the following range: t e {10 6 ; 10 10 }s. In order 
to specify how efficient is the CDI it is sensible to examine an 
accretion process, estimating the corresponding evolution time 
scale, and compare it with that of the CDI. 

Considering the problem of fuelling of AGNs, in 
( King & Pringle 2QQ7\ it has been shown that the self gravitat- 
ing mass in accretion flows can be estimated by the following 
expression: 



M se = 2.76 x 10 5 (—) ^ (—) 
8 \ 0.03/ \0.1) 



-5/27 



0.1 u 



5/27 



M, 



BH 



23/27 



(25) 



where 77 and e are Shakura, Sunyaev viscosity parameter 
( |Shakura & Sunyaev 1973| l and the accretion parameter respec- 




Fig. 6. The behaviour of logarithm of the evolution time scale 
versus M9 and L/Le- The set of parameters is: 77 = 0.03, e — 0.1. 



tively. The latter can be defined by means of the accretion mass 
rate M and the luminosity: 



Mc 2 



Then, defining the accretion time scale t e , 
reduce it to ( |King & Pringle 2007| l: 



(26) 



= M sg /M one can 



'---^"(^kT" 1 ^) 



x 22/27 / £ 

\0.1L £ 



-22/27 



/ M BH 
\ 10 8 M C 



-4/27 



(27) 



As it is clear from Eq. d27l . the accretion evolution time scale 
depends on two major AGN parameters: on the luminosity (L) 
and the AGN mass (Mbh)- Bearing in mind Eq. (|27T i it is sensi- 
ble to investigate t evo i versus L and Mbh- For this purpose let us 
consider the possible maximum "area" of the parametric space, 
L - Mbh studying the behaviour of t evo i for the typical accretion 
disk parameters 77 = 0.03, e = 0.1. 

In Fig.|6]we show the two dimensional surface of logarithm 
of the evolution time scale. The variables are in the following 
range: M g = M B// /(10 9 x M Q ) e {0.01; 1} and L/L E e {0.01; 1}. 
For plotting the figure we took into account that according 
to the observations AGN masses vary in the following range: 
M BH g {10 6 ; 10 9 } x M Q dNelson 20001 . As it is clear from Fig.|6] 
t evo i is a continuously decreasing function of M9 and L/Le- The 
minimum value of the evolution time scale (~ 10 12 s), when the 
accretion process is extremely efficient corresponds to M9 = 1 
and L/Le = 1, whereas the maximum value being of the or- 
der ~ 10 j corresponds to the following pair of the variables: 
M 9 = 0.001 and L/L E = 0.01. 

For understanding how efficient is the curvature drift insta- 
bility, one has to compare the corresponding time scale with the 
evolution time scale. As it has been found, depending on phys- 



ically reasonable parameters t varies in the range: 



10' 



6-10 



whereas the sensible area of t ew i is: ~ 10 12 s. Therefore the 
instability time scale is less than the evolution time scale of the 
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Fig. 7. The dependence of L m /L max versus 6. The set of param- 
eters is: j e0 = {10 5 ; 10 6 ; 10 7 ; 10\R B ~ Ric, n e o = 0.001cm- 3 , 
As = lQQRic, A x = R lc and L = 10 45 erg/s. 



accretion by may orders of magnitude, which means that the lin- 
ear stage of the CDI is extremely efficient. 

The twisting process of magnetic field lines requires a cer- 
tain amount of energy and it is natural to study also the energy 
budget of this process. For this reason one has to introduce the 
maximum of the possible luminosity L max = Mc 2 and compare 
it with the "luminosity" corresponding to the reconstruction of 
the magnetic field configuration L m = AE„,/At w AE„,/t. (here 
AE m is the variation of the magnetic energy due to the curvature 
drift instability). 

Let us consider the AGN with the luminosity, L = 10 45 erg I s, 
then by applying Eq. d26l i and taking into account e = 0.1 one 
can show that the accretion may provide the following maximum 
value: 



10 46 erg/s. 



(28) 



On the other hand if the process of sweepback is realistic, the 
magnetic "luminosity" can not exceed L max . Since this process 
takes place in the nearby zone of the LCS, the effective spatial 
volume can be estimated by the following expression: 



AV * Rf c AR 



4*. 



(29) 



where k = AR/Ri c « 1 represents the non dimensional thick- 
ness of the considered spatial zone. Taking into account Eq. d29l 
the expression of the magnetic "luminosity" reduces to: 



B 2 ^ 

where B r due to the CDI behaves in time as: 
B r = B° r e ,/T 



(30) 



(31) 



Here Z?" is the initial perturbation of magnetic field's toroidal 
component. 

Let us introduce the initial non dimensional perturbation, 6, 
defined as: 6 = B®./Bq, where by Z?q we denote the induction of 
the magnetic field in the leading state [see Eq. (fJTJ]. By consid- 
ering the following set of parameters: y e o = {10 5 ; 10 6 ; 10 7 ; 10 8 }, 



R B * R k , n e0 = O.OOlcm^ 3 , As = l00R k , A x = R, c and 
L = 10 45 erg /s, one can plot the behaviour of L m /L max versus 
the initial perturbation for the characteristic time scale (f w r). 
As we see from Fig. [7J L„,/L max varies from ~ (6 = 0) to 
~ 2.3 x 10~ 9 {6 = 0.1, jeo = 10 8 ). Therefore the maximum 
luminosity (thus the total luminosity budget) exceeds by many 
orders of magnitude the magnetic "luminosity" required for the 
twisting of the field lines. This means that only a tiny fraction 
of the total energy goes to the sweepback, making this process 
feasible. 



4. Summary 

1. Considering the relativistic two component plasma for AGN 
winds the centrifugaliy driven curvature drift instability has 
been studied. 

2. Taking into account a quasi single approach for the particle 
dynamics, we linearized the Euler, continuity and induction 
equations. The dispersion relation characterizing the para- 
metric instability of the toroidal component of the magnetic 
field has been derived. 

3. Considering the proper frequency of the curvature drift 
modes, the corresponding expression of the instability incre- 
ment has been obtained for the light cylinder region. 

4. Efficiency of the CDI has been investigated by adopting four 
physical parameters, namely: the wave length, the flow den- 
sity and the Lorentz factors of electrons and the luminosity 
ofAGNs. 

5. By considering the evolution process of the accretion, the 
corresponding time scale has been estimated for a physi- 
cally reasonable area in the parametric space L - Mbh- It 
was shown that the instability time scale was less by many 
orders of magnitude than the evolution time scale, indicating 
extremely high efficiency of the CDI. 

6. Examining the instability from the point of view of the en- 
ergy budget, we have seen that the sweepback of the mag- 
netic field lines requires only a small fraction of the total 
energy, which means that the CDI is a realistic process. 

An important restriction of the present work is that for de- 
scribing the plasma kinematics a single particle approach has 
been applied. On the other hand it is natural to suppose that the 
collective phenomena must strongly influence the overall kine- 
matic picture of the plasma motion. For understanding how the 
mentioned fact changes the instability, one has to generalize the 
present model for a more realistic astrophysical scenario. 

The next limitation is related to the fact that the magnetic 
field lines were supposed to be quasi rectilinear, whereas in real 
astrophysical situations the field lines might be initially curved. 
This particular case also needs to be studied generalizing the 
present model. 

In this paper the field lines located in the equatorial plane 
have been considered, although in realistic situations the mag- 
netic field lines also might be inclined with respect to the equato- 
rial plane. Therefore it is essential to examine this particular case 
as well and see how the efficiency of the CDI changes, when the 
mentioned inclination angle is taken into account. 
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